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Extrinsic levels, formation energies, and relaxation geome- 
tries are calculated ab initio for oxygen vacancies in a-quartz 
Si02- The vacancy is found to be thermodynamically stable 
in the charge states Q=+3, Q=0, Q =-2, and Q=— 3. The 
charged states are stabilized by large and asymmetric dis- 
tortions near the vacancy site. Concurrently, Franck-Condon 
shifts for absorption and recombination related to these states 
are found to be strongly asymmetric. In undoped quartz, the 
ground state of the vacancy is the neutral charge state, while 
for moderate p-type and n-type doping, the +3 and -3 states 
are favored, respectively, over a wide Fermi level window. Op- 
tical transitions related to the vacancy are predicted at around 
3 eV and 6.5 eV (absorption) and 2.5 to 3.0 eV (emission), 
depending on the charge state of the ground state. 



I. INTRODUCTION 

In this paper we report ab initio calculations of equi- 
librium geometries, formation energies, and extrinsic gap 
levels of vacancies in a-quartz. We concentrate on oxy- 
gen vacancies, which are common defects in both silica 
(amorphous SiC^) and quartz, and alleged to be precur- 
sors of other more complex defects jjj. In view of the 
short-range order conservation upon thermal amorphiza- 
tion of quartz, which basically rules out complex multi- 
center native defects, the state-of-the-art results for oxy- 
gen vacancies in crystalline quartz to be discussed below 
are expected to shed some light on defect phenomenology 
in amorphous silica. Recent ab initio molecular dynamics 
simulations of the amorphization of a small (128-atom) 
sample of a-quartz via melt quenching Q have given ev- 
idence for an almost perfect conservation of short range 
order in the amorphous phase; the structure obtained 
in the simulation Q maintains chemical order, with Si- 
centered tetrahedra linked by corner-sharing oxygens. In 
these simulations, disorder manifests itself merely as a 
broadening of the distribution of Si-O-Si angle between 
linked tetrahedra; the O-Si-0 angle distribution has in- 
stead a sharp maximum close to the tetrahedral angle of 
109°. Edge-sharing tetrahedra are absent in the amor- 
hous phase, but they occasionally appear in the liquid 
On the other hand, even in the liquid the dominant 
defects are simple miscoordinated atoms (3- or 5-fold co- 
ordinated Si, and 1- or 2-fold coordinated O), occurring 
with small (~ 10 %) probability. Although the predic- 
tions of ab initio molecular dynamics may be considered 
too optimistic (mostly, because of the smallness of the 
simulated amorphous sample), it is clear that even the 
simplest defects are not easily formed thermally in SiC>2. 



More complex defects are likely to be even more rare. 
A study of lattice vacancies seems therefore timely, also 
because, although many semiempirical studies |3J have 
appeared in the past, even these simple defects have not 
been widely studied at the ab initio level (recent density- 
functional-theory studies are reported in Q] and Q). In 
particular, no systematic study exists, to our knowledge, 
dealing with the charge states of the O vacancy beyond 
the neutral and +1 states plpf. 



II. METHOD 

We perform ab initio density-functional theory (DFT) 
Q plane-waves pseudopotential calculations, using 
ultrasoft-pseudopotentials @ , which guarantee a high de- 
gree of transferability and a fast convergence in Fourier 
space of all the relevant quantities JtJ. A plane- waves 
basis is used with a cut-off at 20 Ryd. A conjugate- 
gradient total-energy minimization || is used to obtain 
the self-consistent electronic ground-state at fixed ions, 
and an efficient relaxation method j^] is used to reach the 
minimum energy structure following Hcllmann-Feynman 
forces. We employ 36-atom orthorombic supercells with 
no imposed symmetry, and one special k point as Bril- 
louin zone sampling. This setting gives structural prop- 
erties in agreement with those obtained by Liu et al. [[L0| , 
matching experiment to better than 1 %, apart from two 
O internal parameters having a deviation of ~ 2%. 

The formation energy for an O vacancy in SiC>2 in 
charge state Q (where, conventionally, Q is the number 
of electrons transferred from the defect to a reservoir of 
chemical potential /i e ) is given by 



E f {Q) = AE tot (Q) - Mo + Qi^e + E$), 



(1) 



where AE tot (Q) is the total energy difference between 
the defected supercell in charge state Q and the perfect 
bulk supercell, and /i e is the electron chemical potential 
(i.e the Fermi level in our T=0 calculation). The zero of 
He may be assumed to be the top valence band energy 
of the defected crystal in the charge state Q. The 
chemical potential fio is one half of the total energy of 
the spin-unpolarizcd O2 molecule (calculated in a large 
supercell). As expressed above, the formation energy de- 
pends on the charge state and on the Fermi level. The 
latter embodies the effects of background doping, com- 
peting defect species, or externally injected carriers. 
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A. Defect levels as total energy differences 

The eigenvalues of the hamiltonian of a crystal with a 
fixed number of electrons have no rigorous relationship 
to observable optical transition, unless they are many- 
body quasiparticle energies. This is also the case for 
defected crystals, and for transitions involving extrinsic 
levels therein. In particular, transitions to and from ex- 
trinsic levels cannot simply be extracted as differences of 
DFT eigenvalues of e.g. the neutral state of the defect. 
This is true for any one-body theory, including Hartree- 
Fock calculations [jll), since Koopman's theorem gives 
(e.g.) an ionization potential which is not exact even in 
principle. Within DFT one can resort to total energy 
differences, since the DFT total energy is exact in prin- 
ciple, apart from the local density approximation (LP A) 
and other technical issues of practical relevance Q. In 
particular, in a DFT calculation one calculates formation 
(i.e. total) energy differences between the relevant sys- 
tems. For e.g. the case of Q=0 and Q—-1, this difference 
is, by definition, the energy needed to extract an electron 
from the Q=-l state (ionization of the charged center) 
leaving behind a Q=0 center, or to add one to the Q=0 
state (binding to the neutral center) transforming it into 
a Q=-l state. When account for geometrical changes 
upon occupation is taken, this difference is the thermal 
excitation level; "thermal" indicates that the structure 
is allowed sufficient time to rearrange to its new equilib- 
rium position. We concentrate on the calculation of this 
quantity, and discuss later optical transitions. 

Apart from basic choices such as the use of the LDA as 
opposed to gradient-corrected functionals, the uncertain- 
ties in the predicted defect energy levels arise mainly from 
the alignment of the energy zero of differently-charged su- 
percells. Since we choose conventionally the Fermi level 
zero at the valence band top, Q ^ charge states have 
a constant energy offset ~QE® away from the neutral 
state energy (Eq. [l]). We proceed by calculating the bulk 
valence band top E v in a bulk supercell, and express the 
valence top in the charged supercell as 

E$ = E v + AV, (2) 

where AT^ is the difference of the macroscopic averages 
|l3f of the local potentials of the charged supercell and 
of the bulk supercell, 

AV = V Q - F bulk . (3) 

This quantity is evaluated in an appropriate bulk-like re- 
gion, i.e. regions of the simulation cell where the local 
potential is unaffected by the defect. The resulting un- 
certainties are at most of order ±0.1 eV in the thermal 
ionization levels. 



III. RESULTS 

For the neutral state of the vacancy, the formation 
energy is 6.97 eV. Such large energy is in accord with 



the need for irradiation to produce vacancies in the first 
place. As seen from Fig. [y, the vacancy is stable in the 
Q=+3 state for Fermi level values of up to 2 eV above the 
valence band top. Thereafter, the neutral state becomes 
stable, and remains such for Fermi levels up to about 
5.5 eV. The neutral state is therefore the ground state 
of the oxygen vacancy in nominally undoped quartz, i.e. 
for the Fermi level at mid-gap (~ 4.5 eV). Increasing 
further the Fermi level, the vacancy becomes negatively 
charged. The Q — —3 charge state is the most stable 
state for Fermi levels ranging from about 5.5 eV to 8 eV 
(neglecting the small Fermi level range from 5.5 eV to 
5.7 eV, where the -2 state prevails). 

While the thermal level sequence which is usually ex- 
pected is Q = +3, +2, +1, 0, -1, -2, -3 in the present 
case we find a negative-U behavior Ej] between the triply 
positive state and the neutral state, and between the neu- 
tral and the Q = —3 state: the Q =+2 and Q —+1 states 
are unstable against the release of a further electron, and 
the Q=-l and Q=-2 states are unstable towards the cap- 
ture of a further electron. This is due to a strong electron- 
lattice coupling, as is discussed below. The realization of 
large-Q negative and positive states is not unusual for 
wide gap materials. 

A. Defect electronic structure 

The Q = vacancy is a symmetric combination of 
Si dangling bonds. The Si neighboring to the vacancy 
relax almost symmetrically, by 9.5 % of the unrelaxed 
Si-0 bond length, towards the vacant site. In Fig. ||, 
we display density isosurfaces for the defect state of the 
neutral vacancy. 

The only stable positive charge state is the Q=+3 
state. The two Si atoms neighboring to the vacancy 
relax strongly away from their original positions in an 
asymmetric fashion: their distance to the vacancy site 
increases from 1.62 A to 2.16 A, and from 1.61 A to 
2.81 A, respectively. The two Si atom are three-fold co- 
ordinated and close to co-planar with the neighboring 
oxygens. The charge density of the last occupied state is 
delocalized, suggesting that the defect level is resonant 
with the valence band. The thermodynamical instability 
of the singly positive vacancy signals potential difficul- 
ties in previously proposed models Q of the so called E' 
center. 

The Q = — 1 charge state is the antibonding combi- 
nation of the same dangling bonds as above, as can be 
seen in Fig. || (compare with the previous Figure). Ad- 
ditional relaxations with respect to the neutral state are 
below 1%. Thus, electronic promotion into this state will 
require a (rather large) energy of almost purely electronic 
origin. 

In the Q = — 2 charge state, the first antibonding state 
gets completely filled. In contrast to the Q = —1, a 
large, markedly asymmetric relaxation occurs. This can 
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be seen in the structure model of Fig. [|. In this geome- 
try atom Sii has moved towards the vacancy by about 20 
% of the pristine unrelaxed bond length, whereas Si2 has 
moved away from the vacancy by about 3 %. The tetra- 
hedron centered on Sii remains essentially undistorted, 
while that centered on Si2 is distorded concurrently with 
the neighboring tetrahedron centered on a third Si atom, 
S13. This effect is due to a large displacement of the 
bridging O atom. Due to the large relaxation energy gain 
and to the associated electronic rehybridization (visible 
in Fig. |^, where the relevant charge density isosurfaces 
are displayed), the Q = — 2 state is lower in energy than 
the Q = — 1 (Fig]!]) over most of the Ep range, despite 
the electron-electron repulsion. Thus a negative-U effect 
takes place. The -2 state becomes favored over the neu- 
tral at a Fermi level of 5.5 eV. 

At a Fermi level about 0.2 eV higher the Q = —3 
state charge becomes favored. Sii moves towards the va- 
cancy site by an additional 7 % (its distance from the 
vacancy is about 30 % shorther than the unrelaxed bond 
length), while Si2 moves closer by 2 % (its distance from 
the vacancy is about 1 % longer than the unrelaxed bond 
length). The stability of the Q = — 3 state over the 
Q = — 2 state for a wide range of Ep is somewhat puz- 
zling in view of the large repulsion to be expected among 
electrons in the center. The prominent feature of this 
state is indeed that the additional electron minimizes the 
Coulomb repulsion with the preextant ones by occupying 
a state in the interstice between the vacancy-adjacent Sii 
and Si2 atoms, and the S13 atom belonging to the tetra- 
hedron neighboring to Si2 and closest to the vacancy. 
The distance between S12 and Si3 is reduced to only 4.3 
a.u., to be compared with a distance of 5.8 a.u. in the 
unrelaxed structure. This is due to the aforementioned 
distortion of the Si2-0-Si3 bridge linking the two tetra- 
hedral units, whereby the corresponding angle decreases 
from the usual 144° to only about 80°. While this distor- 
tion is in fact already qualitatively present in the Q = — 2 
charge state, the additional electron in the Q = —3 state 
is effectively bound by Si3, as can be seen from Fig. |6|. In 
view of its electronic structure, the -3 charge state may 
be a candidate E'-like center (see below). 

The Q — — 3 state is stable up to extreme n-type con- 
ditions (Fermi level ~ 8.1 eV), whereby the Q = —4 
state sets in. The fourth electron fills in the state whose 
occupation was initiated with Q=-3, with marginal ad- 
ditional relaxation. The charge density and geometry are 
very close to the Q=-3. 

By the time one gets to extreme p or n-type conditions, 
the formation energies of the vacancy become very small. 
Vacancy-driven Compensation of acceptor and donor im- 
purities will be very efficient, and we expect that n-type 
doping of quartz will only be achieved with extreme dif- 
ficulty, if at all. 



B. Franck-Condon shifts 

The Franck-Condon (FC) shift ^Fq,q' is defined as the 
total energy difference between the system in charge state 
Q in its own equilibrium configuration and the same sys- 
tem in the equilibrium configuration of charge state Q 1 : 



(4) 



The FC shift is symmetric upon exchange of Q and Q' 
if the dependence of the total energy on configurational 
coordinates is quadratic and with identical force con- 
stants in the two charge states. However, given the 
large distortions involved in the present case, the FC 
shift might be asymmetric (expecially for large shifts) 
both because of anharmonicity and of force-constant dif- 
ferences. These expectations are indeed verified in our 
calculations. We find large asymmetries for large shifts, 
which are consistent with force constants differing by fac- 
tors of up to about \/2 assuming a single quadratic con- 
figurational cooordinate. We express for convenience the 
FC shifts between the non-positive (non-negative) charge 
states as a matrix jP^eg-states ^os-Btates) whoge ele _ 

ments .^-states ^pos-states-j &re ^ FC ghifts for ^ 

Q = —i (Q = i) state in the geometry of the Q = —j 
(Q = j) state; this matrix would be symmetric in the 
symmetric-FC hypotesis, and its diagonal elements are 
null. For negative states, with Q between and -3, we 
have 



/ 0.60 5.20 6.38 * 
0.37 1.21 1.93 
2.51 1.81 0.18 

\3.17 2.20 0.18 



while for positively charged states (Q between and +3) 
we have 



( 1.35 6.23 
0.64 1.76 
3.78 3.81 0.63 

\4.41 2.90 0.63 



all energies being expressed in eV. Note that only the 
Ti i±i elements are needed in the discussion to follow. 



C. Optical transitions 

The optical transition energy is the total energy dif- 
ference of the two (differently charged) initial and final 
systems in a fixed geometry. Alternatively (but equiva- 
lently), optical levels can be expressed as combinations 
of thermal levels and FC shifts. We only consider one- 
electron (i.e. first-order) processes. This kind of process 
amounts to transferring only one electron at a time into, 
or out of, a defect level. 

The transitions we consider are absorption (an electron 
is transferred from the valence states reservoir into the 
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defect, or transferred from the defect to the conduction 
band via an internal excitation of e.g. the -3 state into 
the -2 state) , recombination (the electron drops into the 
defect level from a reservoir of e.g. optically-excited con- 
duction electrons), and single-particle deexcitation of a 
metastable charged center by recombination of an elec- 
tron into the valence band. We assume that all processes 
involved, e.g. successive absorption transitions (0 — > — 1, 
— 1 — > —2, —2 — > —3 ...), take place on time scales shorter 
than those of lattice vibrations, so that (e.g.) absorption 
occurs in a frozen structure. This assumption seems plau- 
sible in view of typical phonon frequencies and optical 
transition times, especially in large excitation densities. 
We therefore assume that the structure remain frozen in 
the equilibrium configuration of a given state Qf (typi- 
cally the initial state) throughout the transition cascade. 

If the ground state of the defect is charge state Q, in 
terms of thermal levels e[Q/Q'] and FC shifts Fq,q' the 
optical absorption energy for electronic transitions from 
the valence band into charge state Q — 1 is 

<% bs = e[Q/Q - 1] + F Q -i, Qf - F Q , Ql ■ (5) 

For internal absorption excitations (transforming e.g. 
Q = — 2 into Q = — 1 by electron excitation in the con- 
duction band) we have instead 

^t-abs = ^ _ e[Q _ 1/Q] + TQ Qf _ TQ _ 1Qf . ( 6) 

Within the same assumptions, for the conduction band- 
to-acceptor recombination we have emissions at 

^-^gap-^Q 8 (7) 

with -Egap = 8.9 eV, the experimental gap of a— Si02- 
Once again we implicitly assumed that subsequent re- 
combinations are instantaneous on the scale of ionic re- 
laxation. 

Deexcitation is slightly more complicated. Suppose 
that the vacancies have been occupied with a number of 
electrons e.g. between 1 and 4. Assume further that the 
lifetimes of each charge states are larger than the inverse 
of typical phonon frequencies (~ 10 THz). This assump- 
tion seems again admissible in view of typical decay con- 
stants of observed luminescence bands in silica ]l6|] being 
of order 10 -5 to 10 -10 sec. Then, each charge state can 
relax to its own ground state before deexcitation. For 
each vacancy in a given charge state, a cascade of de- 
excitation processes occurs, whereby each charge state 
(between the initial one, say Q — 3, to Q — 0) is occu- 
pied in sequence, and it loses one electron to the valence 
band while in its own equilibrium geometry. Expressing 
again total-energy differences in terms of thermal levels 
and FC shifts, we obtain the dexcitation emission ener- 
gies from charge state Q as 

8%™ = e[Q/Q-l]-T Q , Q -x. (8) 

For internal deexcitation the formula reads 



^nt-dccx - £gap - AQIQ - 1] + FQ-1,Q- (9) 

Here we did not consider selection rules possibly prevent- 
ing some of the transitions involved, and remain with an 
open question on this issue. 

As mentioned above, quantitative predictions depend 
on which ground state is actually realized. This makes 
the optical properties of this system quite rich and com- 
plex. In Table I we summarize the main groups of absorp- 
tion and emission transitions involving different ground 
states. If the Fermi level is at mid-gap (i.e. the solid is 
undoped), the ground state is the neutral state. In mod- 
erate n-type conditions, the ground state will typically 
be the -3 state. The relevant thermal levels are e[0/ — 1] 
= 6.24 eV, e[-l/ - 2] = 5.47 eV, and e[-2/ - 3] = 5.82 
cV. The FC shifts have been given previously. In p-type 
conditions, the ground state will be the +3 state. The 
relevant thermal levels are e[+3/ + 2] = 2.48 eV, e[+2/+l] 
= 2.69 eV, and e[+l/0] = 0.72 eV (the FC shifts have 
been given previously). 

For the neutral state, the absorption band is centered 
around 6.6 eV and the emission band at about 2.3 eV. 
For the -2 charge state there is an absorption band at 
about 6.3 eV, while the emission one is centered around 
3.0 eV. For the Q = —3 state, the absorption bands are 
located around 3.2 eV, 6.0 eV, and 7.1 eV. These corre- 
spond to internal excitations —3 — > —2, —2 — > — 1, and 
— 1^0 respectively. While calculated line intensities are 
not available, we speculate that the latter components 
may give rise to a band centered around 6.3 eV, close 
to an experimentally observed absorption band (6.2 eV) 
usually assumed to be related to the E' center. 

On the other hand, the emission bands are around 2.7 
eV, close to the 2.7/2.8 eV emission band generally at- 
tributed to the oxygen vacancy, and at 3.2 eV, where 
a band is also observed experimentally in silica. For the 
Q = +3 state, absorptions (either normal or internal) are 
predicted at 2.2 eV, 3.1 eV and 5.0 eV, while emissions 
are at 1.9 eV and 0.9 eV. 



D. The Q=-3 as an E'-like center 

We note at this point that the triply-negative charge 
state of the oxygen vacancy exhibits several of the char- 
acteristic features of the E' center || |i"7||l9[ in a-quartz. 
Both the transition generating the centerTjf the starting 
ground state is the neutral vacancy) and that character- 
istic of the center itself are correctly predicted around 6 
eV. While such transition energies appear also in the pre- 
dicted spectra of other possible ground states, only the 
electronic structure of the Q—-3 state seems compatible 
with hyperfme- field and EPR measurements jl^] ■ These 
experiments indicate the presence of two weak hyperfme 
interactions and a strong one associated with the oxygen 
vacancy. In our -3 state the weak signals can be at- 
tributed to Sii and Sia, while the strong signal comes 
from the unpaired electron prevailingly bound to S13. 
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Such attributions could not be made for the Q=+3, 0, 
and -2 ground states. The unpaired electron lobe is ap- 
proximately directed along the direction of the "short 
bond" [p[[l7 in the ideal structure, and as can be seen 
from Fig.TI, S13 is close to three fold coordinated, as it 
is generally accepted that E' should be jl8| . 

The main objection to identifying the Q=— 3 state of 
the vacancy with E' comes from experiments showing 
that E' is positively charged in MOS Si-SiC"2 junctions 
[p7ot| . However, this result is not always reproduced in 
other instances ^0j. It would be interesting to analyze 
further these experiments in the light of our results, tak- 
ing into account for instance band offset effects at the 
Si/SiOo interface. As to magnetic data, recent calcula- 
tions |5[| of the hyperfine tensor of a distorted variant of 
the singly positive vacancy have resulted in good agree- 
ment with the data pertaining to the E' center JltJ , giving 
support to the attribution of E' to a positive charge state 
of the vacancy. On the other hand our results indicate 
that the Q=+l vacancy is unstable towards becoming 
triply charged, and that the Q=+3 state appears to be a 
valence resonance, lacking features which could explain 
the peculiar magnetic signature of E'. 

It is clear that additional work is needed to resolve 
this puzzle. Nevertheless, the 3- oxygen vacancy seems 
to carry several E'-like features, and it may be considered 
as a possible new variety of E' defect. 



IV. SUMMARY 

We have presented ab initio results on oxygen vacancies 
in a— SiC>2. The ground states for the vacancy as func- 
tion of background doping are the triply-charged positive, 
neutral, and triply-charged negative states. The high- 
charge states are stabilized by large, asymmetric distor- 
tions involving at least three tetrahedral units neighbor- 
ing to the vacancy. In connection with the possible differ- 
ent ground states we predict optical absorption energies 
in the 3.2 eV and 6-7 eV ranges, and emissions in the 
2.3-3.2 eV range. Based on predicted optical transitions, 
geometry, and electronic structure, we suggested a possi- 
ble new E'-like center, namely the triply-charged negative 
oxygen vacancy. 
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Ground state 


absorptions 


emissions 


Q = +3 
Q = 
Q = -2 
Q = -3 


2.2 / 3.1 / 5.0 ± 0.2 

6.6 ± 0.2 

6.3 ± 0.3 
3.2 / 6.0 / 7.1 ± 0.2 


1.9 / 0.9 ± 0.1 
2.3 ± 0.2 
3.0 ± 0.1 
2.7 ± 0.2 / 3.2 ± 0.2 



TABLE I. Summary of optical transitions (in eV) associ- 
ated to the oxygen vacancy in a— Si02 
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FIG. 1. Formation energy of the oxygen vacancy in a-Si02 
as a function of the Fermi level. 

FIG. 2. Charge density isosurfaces for the neutral state 
of the oxygen vacancy. The surface shown corresponds to a 
density value of 0.02 electrons/bohr 3 . 

FIG. 3. Charge density isosurfaces for the singly-negative 
charge state of the oxygen vacancy. The surface shown corre- 
sponds to a density value of 0.006 electrons/bohr 3 . 

FIG. 4. Relaxed geometry of the Q = —2 state of the oxy- 
gen vacancy (black atoms and thick connecting lines). The 
geometry of the neutral vacancy is given for reference (white 
atoms and thin connecting lines) . 

FIG. 5. Charge density isosurfaces for the doubly-negative 
charge state of the oxygen vacancy. The surface shown corre- 
sponds to a density value of 0.02 electrons/bohr 3 . 

FIG. 6. Charge density isosurfaces for the triply-negative 
charge state of the oxygen vacancy. The surface shown corre- 
sponds to a density value of 0.0045 electrons/bohr 3 . 
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